Interaction induced phase fluctuations in a guided atom laser 
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In this paper, we determine the magnitude of phase fluctuations caused by atom-atom interaction 
in a one-dimensional beam of bosonic atoms. We imagine that the beam is created with a large 
coherence length, and that interactions only act in a specific section of the beam, where the atomic 
density is high enough to validate a Bogoliubov treatment. The magnitude and coherence length of 
the ensuing phase fluctuations in the beam after the interaction zone are determined. 



Recent progress on micro-fabricated magnetic traps 
and guides for atoms |) have stimulated the efforts 
to realize guided atom-laser beams || with possible ap- 
plications for high precision atom interferometry. The 
phase coherence is a crucial property of such beams and 
interferometers, and interactions between atoms may in- 
fluence the coherence. So far, the effect of interactions on 
the condensate phase were mainly studied in condensates 
of finite spatial extent Q ,where both degradation || and 
squeezing of the phase g have been studied. In these 
cited works, the main effect can be ascribed to the num- 
ber fluctuations in the condensate and the dependence of 
the mean field energy on atom number due to interac- 
tions. The goal of the present work is to determine the 
amplitude and the coherence length of the phase fluctu- 
ations created by the atomic interactions in a continuous 
beam whose coherence length is initially infinite. The 
main contribution in this case is collisions that transfer 
atoms from the condensate at momentum K to different 
momenta K + k in the beam. At very low density an 
interacting bose gas in one dimension is described by the 
Tonks-Girardeau regime, but we shall assume a higher 
density for which Bogoliubov theory is well suited to de- 
scribe the deviations from mean field theory @. 

We consider the situation depicted in Fig.l, where a 
non-interacting mono-energetic beam of atoms with mo- 
mentum TiK arrives in a region of length L where the 
atoms interact with each other. For the practical real- 
ization of such a model, one may imagine a guided atom 
beam, where the transverse width is particular narrow, 
and therefore the density and the effective interaction is 
significantly increased over the length L. 

The system is described by the effective one dimen- 
sional Hamiltonian 

H = j dxif/t (x)hQip(x) + |j J dxip^ (x)^ (x)ip(x)ip(x), 

(1) 

where ho = —h~/(2m)d 2 /dx 2 , and the coupling constant 
g = ^-|r, assuming that the size of the transverse 

ground state in the guide a± = y /, 2h/muj is much larger 
than the s-wave scattering length a. 

We always consider the case where na\/a >> 1 so that 
correlations between atoms are small, and the state of the 



FIG. 1: An atomic beam crosses a region L where the atoms 
interact with each other. The interactions scatter atoms out 
of the macroscopically populated mode into other momentum 
states, and the resulting phase fluctuations in the beam can 
be probed by the interference with a coherent reference beam. 



system can be described by a mean field wave function 
and a noise term 



ip(x) 



iKx 



5ip(x). 



(2) 



The fluctuations represented by Sip(x) are probed by 
interference with a reference beam in a coherent state 
ip re f(x) = ae lKx . The quantity measured is the num- 
ber of atoms on each side of the beam splitter. In an 
actual experiment, this beam splitter may be realized by 
tunneling between two guides, or the atomic beams may 
be physically overlapping, but with two different internal 
states, which can be coupled by Raman laser beams. In 
either case, the operators giving the density of atoms in 
the two output states are 

n±(x) = \{i>l ef {x) ± tf{x)) (Aef(x) ± tl>(x)) (3) 

For a = Y / n u "e I7r / 2 , the mean density is identical in the 
two output beams and we define the local phase operator 
as 



0(x) 



n + (x) — 7i- (x) 



(4) 
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The fluctuations in Q(x) are given by the correlation func- 
tion 

{9{x)9{x')) = - {-{8xj;(x')5tp(x))e- iK( - x+x ^ - c.c. (5) 
" +{5^{x')5il>{x))e ui ^'-^ +c.c. 
+S(x' - x)}, 

and it is our goal to present a quantitative analysis of 
these fluctuations. After averaging over a distance X, 
the fluctuations of the phase is given by 



Ae -\\xL 6{x)dx 



(0) 



In the absence of interactions between atoms, only the 
last (^-correlated quantum noise term is present in Eq.(|5|), 
and we find AO = \j\fnX as expected from the usual 
number and phase uncertainties in a coherent state. 

In the presence of interactions between atoms, a con- 
densate will experience phase diffusion which is due to 
the spread of chemical potential over the Poisson distri- 
bution of number states which form the coherent state. 
Let us assume that the atomic beam is in fact a very long 
wavepacket with length C. When a superposition of num- 
ber states | TV) of this wavepacket passes the interaction 
zone the interaction perturbs the energy of each number 
state by the amount \n 2 L = &^L. The passage time is 
T = C/(hK/m) and hence the number state component 



| TV) experiences a phase shift On 



grtiL 
2Ch' 2 K 



TV . The phase 



of a coherent state is the derivative dO/dN of the phase 
of the Fock components. After passage of the interacting 
zone, this derivative depends on TV and the spread of the 
phase over the width y/~N of the Fock state distribution 
is 



AO 



gmL ^JngmL 

Vch 2 k' 



Ch 2 K 



(7) 



The relative number fluctuations vanish in the limit of 
very large coherence length £ of the incident beam, and 
in that limit the Bogoliubov excitations will give the dom- 
inant contribution to phase diffusion. 

If the expression for the field operator (Q) is introduced 
in the second quantized Hamiltonian (|l|), and the noise 
terms 6ip(x),8i/j' (x) are truncated above second order, 
H can be diagonalized by a Bogoliubov transformation, 
H = Eq + ^2<jJka\ak, where the field operator 5ip(x) is 
expanded on the bosonic operators oik,oc\ as Sip(x) = 
J2k (uk(x)ak — vl(x)a\^j. The wavefunctions u k and Vk 
solve the equations 



h - fi + 2g\ip\ 2 -gtp 2 \ f u k \ _ ( u k 

gip* 2 -h Q + u-2g\<p\ 2 J [v k J -^{vk 

(8) 

where g takes the value zero outside the interaction in- 
terval, and u k and v k fulfill the normalization condition 

/ dx(u k (x)uk>(x)* - Vk(x)vk> (x)*) = 8(k - k'). 



We neglect reflection of the mean field wave function at 
the entrance and exit of the interaction region, which is 
a good approximation as long as gn <C h 2 K 2 /2m. When 
we write the mean field wave function in the interaction 
zone as ip(x) = s/rie , the chemical potential is \i = 
h 2 K 2 /2m + gn, and Eq.(@) has the plane wave solutions 



e i(k+K)x jj 
e i(k-K)xy 



(9) 



Inserting this expansion in Eq.(^|) we obtain 

- fx + 2gn\ U - gnV = ujII 

- (i + 2gn) V - gnU = -wV 




(10) 



Outside the interaction zone, where g = 0, U and V are 
independent with energies u) u i v = kK ± k 2 /2. Inside the 
interaction zone where g ^ 0, we multiply the first equa- 
tion of Eq.(JlO|) by (j^(k - K) 2 - u + 2grij and obtain 
a second order equation for oj which leads to 



lu = — kK ± 

m 



<h 2 k 2 fn 2 k 2 



2m V 2m 



2gn 



(11) 



The spectrum has two branches corresponding for high 
energies to particle and hole-like excitations, respectively. 
For the two branches, the ratios V/U of the solutions are 



2/2 



ii 



in t/n 

h 2 

mgn 



k 2 , k' 2 
ZLl —L 
2 V 2 



2mgn/h 2 ) 



L (%- + 2mgn/h 2 ) 



(12) 

In each region along the atomic beam, four solutions 
are possible, corresponding in the noninteracting region 
to particles or holes propagating in one direction or the 
other, and the global solutions should be continuous with 
continuous derivatives for both the u and v functions. 
Thus 8 equations relate the 12 parameters and for each 
energy 4 solutions exist. Like in normal scattering theory, 
our choice of boundary conditions serves to identify the 
relevant basis for the solution to the problem. Excitations 
with an escaping hole (removal of an incident particle) in 
either direction are not physically relevant, and we thus 
study the elementary excitations with incoming particle 
components only. 

Let us first consider the excitation corresponding to an 
incoming particle u(x) oc e l ( K+k ) x from the right, (k < 
—K). This particle has a momentum which differs from 
that of the condensate by more than K . If we assume that 
rather than being an exact delta-function, the atomic in- 
teraction has a momentum cut off which is smaller than 
K, the coupling to the v(x) function in the Bogoliubov 
equations is suppressed, and this mode does not con- 
tribute to phase fluctuations in the Bogoliubov vacuum. 



3 



Another set of Bogoliubov modes corresponds to incom- 
ing particles from the left with momentum k + K and 
amplitude Ui — -y= (the condition [5ip(x), 8ip'(x')] = 
8{x — x') serves to normalize the Bogolibuv mode since 
J k dkU l {k)e^ K+ ^ x U. l {k)*e^ K+k ^' = 5(x - x') on the 
left side where the Bogoliubov modes contain no holes). 
Eq.(|lT|) has four solutions, but wc will consider only mo- 
mentum components K + ki and K+ fe close to K inside 
the interaction zone, since the effective momentum cut-off 
in the interaction prevents coupling to waves with very 
different momenta. Continuity of the solution at the left 
entrance to the interaction zone implies 



U i = U 1 + u 2 
= Vi + v 2 



(13) 



Eq.(|l2|) dictates the ratio between the U and V ampli- 
tudes, and the continuity of the derivatives of u and v 
cannot be fulfilled without allowing for reflection. This 
reflection will, however, be very small as k <C K, and we 
make an insignificant error by requiring only continuity 
of the Bogoliubov mode functions. All amplitudes in the 
inteaction zone are hence given by the incident amplitude 



1 l V2-Vi anr | J 1 2 1/2-2/1 (1A\ 



■ VI— 



■ 1/1—2/2 



At the right of the interaction zone, we have to match to 
the independent solutions jj e l ( K + k ) x and Ve~ l ( K ~ k ) a of 
same Bogoliubov energy to, i.e., the amplitudes U and V 
are given by the continuity relations 



Ue i(k+K)L _ V ^ e i(k x + K)L + jj 2e i(k 2 +K)L 
y e i(k'-K)L = y ie i(k 1 -K)L + y^i^-K)! 



(15) 



Using the relations (!#, this gives 



U — JJ. e —ikL I V2 ik-\ L , _m_ e ik 2 L 
1 \V2-yi 1/1-1/2 / 

V = U ie -* k ' L (jU^) e l ^ L 2ism(^L) 

The phase fluctuations given in Eq.(^|) can now be ex- 
pressed as 



(6(x)8(x')) = - dk {2Ke(V(k)*U{k)t 



a —ik 

~2\V(k)\' 2 cos(k , {x / -x))}, 

(17) 

where the (5-function part has been omitted. 

We now present a number of simplifications, which are 
valid under reasonable assumptions and which will lead 
to a better understanding of the dependence of the phase 
fluctuations on the physical parameters of the problem. 
First, as long as ki ~ k 2 — k <C K we can use the 
expressions of y\ and y 2 ( p^| ) for the single value k. We 
cannot, however, replace k\ and k 2 generally by k in the 
exponential factors in (16). Instead, we look at limiting 



cases. If k ^> y/mgn/h, the energy oj = kK + k 2 /2 leads 
to the same dispersion law as that of a non interacting 
gas. This dispersion law leads to the expressions, hi ~ 

,2 

k + ^ and k 2 ~ fc, and the expressions for y\ and y 2 
are in this limit approximated by y\ ~ an d yi 

le V(k) amplitude in Eq.(! 
<C 1 results in only negligible 



WW 

by the factor 



j-Tj^- The suppression of the V(k) amplitude in Eq.(16) 



mj/2 



mgn 



contributions to phase fluctuations from these momentum 
components. If instead, \k\ <C ^/mgn/h, the Bogoliubov 

dispersion law approximately yields fci ~ k + k ^"5" and 



UK 



KK 



and the expressions for yi and y 2 are 



in this limit approximated by y\ 



and y 2 



1 



h\k\ 



We then find, 



U{k) ~ V(k) 



lkL o ' V rn 9 n "'ki+k2 

7 2 ^ Zl ^Kk~ 



(18) 



In Eq.(p7|) k' is the momentum associated with the v- 
componcnt in the nonintcracting region and, as long as 
|fc| <C K, it writes k' ~ k + j^, and the cross term in 
V(k)*U(k) can thus be written, 

2 



I(x,x') 



dk\V{k)\ 2 Ke{e lk{x '- x) e lxk2/K ). (19) 



Approximating k' by k the |U(fc)| 2 term in Eq.(0) can 
be written 

J(x,x') = - [ dk\V(k)\ 2 cos(k(x- x 1 )). (20) 
n J 

For large x, I(x,x') becomes wider as a function of 
y = (x' — x). If wc label x/K as t/2m in the last 
exponential we recover the expression for the spreading 
wave packet of a massive particle, given by the momcn- 



hx 



turn spread Afc of |F(fc)| 2 , i.e., Ay = hAkj^ - 

At distances far from the interacting region I(x, x') gives 
thus a negligible contribution to local phase fluctuations, 
but when integrated over intervals larger than Ay, the 
conserved 'norm' of the spreading wave packet implies a 
result comparable to the contribution from the |U(fc)| 2 
term, estimated below. 

Fig. 2a). shows the value of (6(x)9(x')) for different 
values of £ and as functions of x' . The curves are obtained 
by a numerical integration of Eq. (|l7j) with the proper 
expressions for U(k) and V(k). For large x, we both see 
the spreading I(x 7 x') and the narrow J(x 7 x') component, 
cf. Fig 2b). 

We are now in position to obtain the essential scaling 
of the phase fluctuations with the physical parameters of 
the problem: If wc focus on the narrow part J(x, x'), at 
distances from the interaction zone larger than L 2 gn/K, 
the phase fluctuation have a finite coherence length, esti- 
mated as the inverse of the width in momentum space of 
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V(k), 



Ax 



y/mgnL 
hK ' 



(21) 



The dependence on the size of the interacting region L 
is due to a phase matching conditions : momentum and 
energy conservation cannot be both fullfilcd in a collision 
which creates excitations and the violation of momentum 
conservation becomes more severe as L increases. The 
amplitude of the fluctuations increases linearly with L, 



(9{xf) 



i / \ 3/2 T 

1 / mgn \ L 
n (j? 



K 1 



(22) 



and to recover a well defined phase, one has to integrate 
over a finite detector region in space. If one detects atoms 
on a length X larger than Ax, the phase precision is given 
by the k = component of the fourier transform of the 
phase fluctuations, 



A 2 = ^ |y(o)|2= - 
n nX 



2 / mgnL\ 



\ h 2 K ) 



, W 2 ) 

' X/ Ax 



(23) 



If X becomes larger than the length scale Ay, the integral 
of J(x, x') contributes with the same amount, i.e., Eq.(E3|) 
is multiplied by a factor of 2. This analysis is confirmed 
in Fig. 2c)., showing as a dashed line the analytical result 
( p3| ) and as a solid line the results of the full numerical 
calculation. The dot-dashed line in the figure shows the 
value (22), applicable for short intervals X < Ax. 

Let us conclude with a numerical example for a beam 
of Rubidium atoms with velocity TiK/m = 80mm/s, and 
a linear density of n = I0 6 atoms/m, subject to an 



interacting region with a perpendicular confinement of 
a± = 0.5 /jm and a length of L = 1 cm. With the value 
a = 5 nm for the scattering length we then obtain 



Ax ~ 26 fj,m and (8(x) 2 ) ~ 2. 



(24) 



To have a phase uncertainty AO as small as 0.1, one has 
to count the Rb atoms in this example over a length of 
order 2 cm. Note that the phase uncertainty is two orders 
of magnitude larger than the standard result for a per- 
fect coherent state on that interval, and if we assume a 
coherence length C of the beam exceeding 2 cm,, the con- 
tribution from number fluctuations (7) is smaller than the 
Bogoliubov excitation result by at least a factor of two. 

To summarize, we have quantified the phase fluctu- 
ations in an atomic beam, and found that they mani- 
fest themselves at different observational length scales: 
If X 1/n, the number of atoms observed is of order 
unity and the phase is always ill-defined. The phase un- 
certainty falls off when X is increased, but in the pres- 
ence of interactions, until X reaches the length scale 
Ax = y'mgnL/hK, it levels off to a value given by 
Eq.([22]). For X larger than Ax, it decreases according 
to Eq. (|23|) or Eq.(|23|) multiplied by a factor two if X is 
also larger than Ay = Kxq / (L^Jmgn) , where xo is the lo- 
cation of the detection interval. If the coherence length C 
of the beam is short, there is an extra contribution (j7|) to 
phase fluctuations. Finally, it should be mentioned that 
other sources of phase fluctuations may be important in 
a guided atomic beam. For example, the decoherence 
and heating of the atomic beam due to the presence of 
the surfaces of the macroscopic elements that provide the 
guiding potential has been analysed in ||. 
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FIG. 2: Phase fluctuations calculated for an exanmple with numerical parameters L = 10 4 /K and gn = 0.02h 2 K 2 jm. (a) : Phase 
fluctuations given by Eg. (f7?|). For large x the curves clearly separate in a central feature and a broad pedestal due to the \V(k)\ 2 
and the V(k)*U(k) contributions, respectively, (b) : Phase correlations around xo = 20L 2 gn/ K . The J(x,x') term is shown as 
a dashed line, (c): Phase fluctuations averaged over intervals of variable length X. The curves show dx6{x)/X) 2 ) for 

xo = 50L 2 gn/ K . The full line is the result of a numerical calculation, the dashed line is the analytical expression (fji|j, and the 
horizontal dot-dashed line is the value $2!^ ) expected for X smaller than the coherence length Ax = ^/gnmL/hK . 



